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We use QCD sum rules to study the masses of the baryons H c and Hi,. We work with a current 
where the strange and the light quarks are in a relative spin zero, at leading order in a s . We consider 
the contributions of condensates up to dimension six. For Ha we get ms b = (5.75 ± 0.25) GeV, and 
for H c we get m= c = (2.5 ± 0.2) GeV, both in excelent agreement with the experimental values. We 
also make predictions to the state Qt,(ssb) obtaining mn b — (5.82 ± 0.23) GeV. 
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The recent observation of the baryon by DO and CDF collaborations [l|, l2| with mass in agreement 
with the prediction in ref. [H, has estimulated us to use the QCD sum rules (QCDSR) [J] to evaluate the 
mass of this state. Previous QCD sum rule calculations for beauty baryons have been done before for 
the Ab, Sf,, and for double beauty baryons [j| but not for S^(ds6). Since in the QCDSR approach, 
hadronic masses are related with the vaccum condensates, the use of this method to evaluate hadronic 
masses is an important step in the understanding of the dynamical nature of these masses. 

Here we follow ref. [3J] and assume that the strange and light (sq) quarks in are in a relative spin zero 
state. Therefore, the most general (low dimension) current which interpolates the operator can be 
constructed from a combination between the two currents, formed with scalar and pseudoscalar diquarks: 

VQ = e abc [(qICj 5 s b ) + t(q^Cs b )^ 5 ] Q c , (1) 

where a, 6, c are color indices, C is the charge conjugation matrix, Q denotes the heavy quark and t is 
the mixing parameter between the two currents. Of course the above interpolating field can also be used 
to study the S c (qsc) baryon. 

The QCDSR is constructed from the two-point correlation function 

11(g) =ijd 4 x e^(0\T[r, Q (x)fj Q (0)}\0)- (2) 

Lorentz covariance, parity and time reversal imply that the two-point correlation function in Eq. ([2]) has 
the form 

n{q) = n 1 {q 2 )+0l 2 {q 2 ). (3) 

A sum rule for each invariant function II i and I±2, in Eq.([3]) can be obtained. 

The calculation of the phenomenological side at the hadron level proceeds by writing a dispersion 
relation to each one of the invariant functions in Eq. ([3]): 

nf e "(g 2 ) = -/rf g 2 Pl{s \. +■•-, (4) 

J q z — s + te 

where pi is the spectral density and the dots represent subtraction terms. The spectral density is de- 
scribed, as usual, as a single sharp pole representing the lowest resonance plus a smooth continuum 
representing higher mass states: 

Pi(s) = \ 2 m SQ 5(s-m% Q )+p? nt (s), 

P2(s) = X 2 5(s-m% Q ) + p c 2 ont (s), (5) 
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where A 2 gives the coupling of the current with the low mass hadron of interest. For simplicity, it is 
assumed that the continuum contribution to the spectral density, pf nt (s) in Eq. (O, vanishes bellow a 
certain continuum threshold sq. Above this threshold, it is assumed to be given by the result obtained 
with the OPE. Therefore, one uses the ansatz Q 

P? nt (s) = P ? PE ( s )e(s - s ) , (6) 

with 

p ope [s) = I 7m[n ?^( s )] . (7) 

On the OPE side, we work at leading order in a s and consider the contributions of condensates up to 
dimension six. We keep the terms which are linear in the strange-quark mass m s . We use the momentum 
space expression for the charm quark propagator, while the light-quark part of the correlation function 
is calculated in the coordinate-space. After making a Borel transform of both sides, and transferring the 
continuum contribution to the OPE side, the sum rules for Sq baryon, up to dimension-six condensates 
can be written as: 

f So 

\ 2 m 3Q e~ m ^ 2/M2 = ds e~ s ' M2 p? PE (s) + H 1 (M 2 ) , (8) 

Jm 2 Q 



A e _ mH< 5 



dse-^ 12 p OPE (s) +n 2 (A/ 2 ), 



where 



(S) = ffrr is)+p m [s)+prl{s) 

n,(Af 2 ) = nf 9a - Gq) {M 2 ) + iif q)2 {M 2 ). 



In the structure 1 we get 

(l-t 2 )m 5 r 



2 7 7T 4 

mom. 



(1 - x) ( i + — + 1 ) + 6 ( 1 + ~ ] In a; 



-(1-x) 



2 3 7T 2 

(l-i 2 )m Q (. 9 2 G 2 ) 



(l-t 2 )Q 



2 9 3tt 4 



(1 -x) 7 + - +6 In a; 



(qga.Gq) ( A ,2, _ m Q m s 



n\ qgaKjq > (m 



Iif q)2 (M 2 ) 

where x = ttiq/s. In the structure ^ we get 
(\ + t 2 )m% 



2 5 7 



(l-t 2 )(sga.Gs) e _ ml/M . + t 2 ){ - g(jGq) (-m^M* 



™q(gg)B (1 , t 2 )e -m%/M^ 
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(10) 



(11) 



2 5 7T 2 



(1 + 1 )(sga.Gs) e _ m%/M , + (1 _ t 2 ){ - ga Gq) L-<I m2 + J^ da(l - a)e<^ 



M5(l_t2 )e -m 2 Q /M 2 _ 



(12) 
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The contribution of dimension-six condensates (g 3 G 3 ) is neglected, since it is assumed to be suppressed 
by the loop factor 1/167T 2 . 

In the numerical analysis of the sum rules, the values used for the quark masses and condensates are 
(see e.g. 0, ||): m s = (0.10 ± 0.03) GeV, m c (m c ) = (1.23 ± 0.05) GeV, m b (m b ) = (4.24 ± 0.06) GeV, 
(qq) = -(0.23±0.03) 3 GeV 3 , fi = 0.8(qq), (qga.Gq) = m 2 (qq) with to 2 = .8 GeV , (g 2 G 2 ) = 0.88 GeV 4 . 

We start with the charmed baryon E c . We evaluate the sum rules in the range 1.5 < M 2 < 3.0 GeV 2 
for so in the range: 3.0 < ^/sq < 3.2 GeV. In Fig. [T] we show the contribution of each term in Eq. (fT2|) to 
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FIG. 1: The OPE convergence for the sum rule Eq.© for S c , using y/so = 3.1 GeV and t = 1. The dotted, 
long-dashed, dashed and dot-dashed lines give, respectively, the perturbative, quark condensate, gluon condensate 
and mixed condensate contributions. The solid line gives the total OPE contribution to the sum rule. 

the sum rule in Eq. |9j) , for t = 1 and y/so — 3.1 GeV. We see that we get an excelent OPE convergence. 
For t = 1 the four-quark condensate contribution to the sum rule Eq. (0) vanishes. For other values of t, 
although the four-quark condensate contribution is bigger than the contribution of the other condensates, 
it is still much smaller than the perturbative contribution and, therefore, it does not spoil the convergence 
of the sum rule, as can be seen in Fig. [2l Therefore we conclude that the convergence of the sum rule 
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FIG. 2: Same as Fig. 1 for t = 0. The solid line with dots gives the four-quark condensate contribution 

in the $ structure in Eq. @, is good for any value of the mixing parameter t (see Eq.([T])) in the range 
< t < 1. This is not the case of the sum rule in Eq. © (structure 1), since for t = 1 the perturbative 
and gluon condensate contributions vanishe. As a matter of fact, if we try to obtain the mass of the S c 
baryon by dividing Eq. by Eq. ©, we only get values compatible it the experimental mass for t ~ 0. 
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Therefore, in this work we will use only the sum rule in Eq. To obtain the mass of the baryon we 
take the derivative of Eq. ([9]) with respect to 1/M 2 , divide the result by Eq. ([9]) and obtain: 



ml 



H ds e- s ' M * 



s p° PE (s) 



{dn 2 /dM- 2 ) 



CI ds e 



-s/M 2 p OPE^ 



n 2 (M 2 ) 



(13) 



We get an upper limit constraint for M 2 by imposing the rigorous constraint that the QCD continuum 
contribution should be smaller than the pole contribution. The maximum value of M 2 for which this 
constraint is satisfied depends on the value of sq and t. The comparison between pole and continuum 
contributions for y/so = 3.1 GeV and t = 1 is shown in Fig. [3J The same analysis for the other values of 




M (GeV ) 



FIG. 3: The dashed line shows the relative pole contribution (the pole contribution divided by the total, pole 
plus continuum, contribution) and the solid line shows the relative continuum contribution for y^so = 3.1 GeV 
and t — 1. 



the continuum threshold and t = 1 gives M 2 < 2.3 GeV 2 for y/s~E = 3.0 GeV and M 2 < 2.65 GeV 2 for 



y/so = 3.2 GeV. We get similar results for other values of t, for example for ^/sq = 3.1 GeV and t = 
we get M 2 < 2.6 GeV 2 . 
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FIG. 4: The H c baryon mass as a function of the Borel parameter (M 2 ) for different values of the continuum 
threshold and the current mixing parameter: t = 1, t/so = 3.2 GeV (solid line); t = 1, ^/sio = 3.0 GeV (dotted 
line); t = 0, ^/Iq = 3.1 GeV (dot-dashed line). 



In Fig. 01 we show the S c baryon mass obtained from Eq. (|13p . in the relevant sum rules window for 
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different values of ^/so and t. From Fig. 0] we see that the results are more stable, as a function of M 2 , 
for t = 1 than for t = 0. Therefore, we will use t = 1 to estimate the mass of the particle. It is very 
interesting to notice that the result for the mass obtained with t = is very similar to the one obtained 
dividing Eq. © by Eq. © using also t = 0. 

We found that our results are not very sensitive to the value of the charm quark mass, neither to 
the value of the condensates. The most important source of uncertainty is the value of the continuum 
threshod and the Borel interval. Using the QCD parameters given above, the QCDSR result for the S c 
baryon mass is: 

m Sc = (2.5 ±0.2) GeV, (14) 

in a very good agreement with the experimental value m^ v = (2.4710 ± 0.0004) GeV 0]. 

In the case of the beauty baryon 3j,, using consistently the perturbative MS-mass mb(mf,) = (4.24 ± 
0.6) GeV, t = 1 and the continuum threshold in the range 6.3 < y/so < 6.5 GeV, we find a good OPE 
convergence for M 2 > 4.0 GeV 2 . We also find that the pole contribution is bigger than the continuum 
contribution for M 2 < 5.2 GeV 2 for < 6.3 GeV, and for M 2 < 5.7 GeV 2 for = 6.5 GeV. 
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FIG. 5: The H;, baryon mass as a function of the Borel parameter (M 2 ) for different values of the continuum 
threshold and the current mixing parameter: t = 1, y/so = 6.3 GeV (solid line); t = 1, = 6.5 GeV (dotted 
line); t = 0, ^/so = 6.4 GeV (dot-dashed line). 

From Fig. where we show the 5& baryon mass, we see that the results are very stable as a function 
of M 2 in the allowed Borel region. For completeness we also show, in Fig. [5l the results obtained using 
t = and y/so = 6.4 GeV. In this case the results are very stable for all values of t in the range < t < 1. 
Therefore, we will also use different values of t to estimate the uncertainties in the result. Taking into 
account the variation of M 2 and varying sq, t and mi, in the regions indicated above, we arrive at the 
result: 

m Eb = (5.75 ±0.25) GeV , (15) 

also in a very good agreement with the predictions in refs. Q and Q, and with the experimental results 
in ref. l\: m§° = (5.774 ± 0.013) GeV, and in ref. 0: mg^ F = (5.7929 ± 0.0024) GeV 

We have presented a QCDSR analysis of the two-point functions of the Sg(qsQ) baryons. We find 
that the sum rules results for the masses of S c and are compatible with the experimental values and 
with the predicitons in refs. [1, Q, in the case of S&. These results for Hj are sumarized in Table I. 



Table I: Results for the 5^ mass. 



m Sb (GeV) 


ref. 


5.774 ±.011 ±0.015 


DO [1] 


5.7929 ±0.0024 ±0.0017 


CDF [2] 


5.795 ± 0.005 


[3] 


5.8057 ±0.0081 


[8] 


5.75 ±0.25 


this work 



G 



It is important to notice that while the calculation based on modeling hyperfine interaction [H, Q 
needs, as inputs, the masses of others heavy baryons, in our calculations the masses are extract using 
only information about the QCD parameters as quark masses and condensates. In the case of the heavy 
baryons S c and in particular, their masses are determined basicaly by the first term in the sum rules 
Eqs. (J5J and © and the condensates are not very important. 

We have tested two different choices of currents. While the results for S c are sensitive to this choice, 
in the case of Sb the results are very stable for the mixing parameter in the range < t < 1. 

It is not possible to generalize directly our results to the baryons Q c (ssc) and Qb(ssb), from the current 
in Eq.{T]), since one can not construct an scalar or pseudoscalar diquark in a 3 configuration of color, 
with two strange quarks. Therefore, to study the baryon Qq we can use a proton-like current: 



(16) 



With this current, the OPE contributions to the sum rule in the structure $ (Eq. ([9])) for the £Iq baryon, 
up to dimension-six condensates, are given by: 
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(17) 



We find that the OPE convergence of this sum rule is as good as the OPE convergence obtained for 
E Q in the same Borel region: M 2 > 1.5 GeV 2 for fl c and M 2 > 4.0 GeV 2 for Q b . In Table II we give 
the maximum value of M 2 for which the continuum contribution is smaller than 50%, for different values 
of the continuum threshold. 

Table II: Maximum values of the Borel parameter (M 2 ) for different values of So- 



state 


GeV) 


ML,. (GeV 2 ) 


n c 


3.2 


2.7 


n c 


3.4 


3.1 


n b 


6.4 


5.5 


n b 


6.7 


6.5 



The results obtained for mn Q , from Eq. lfT3")) . are also very stable, as a function of M 2 , in the allowed 
Borel window. Considering the variations in sq and M 2 given in Table II, and the variations in the quark 
masses and condensates, as discussed above, we get 



mn c = (2.65 ±0.25) GeV , 



(18) 



in a very good agreement with the experimental value m^f = (2.6975 ± 0.0026) GeV 0. For fij, we 
make the prediction: 



(5.82 ±0.23) GeV . 



(19) 



As a final remark, it is very reasurring to see that the OPE convergence is so good for heavy baryons, 
since this is not the case for tetraquark states Q. As shown in ref. [§j], it is very difficult to find a 
Borel region where the continuum contribution is bigger than the pole contribution and where the OPE 
convergence is acceptable, for tetraquark states with only one heavy quark. Therefore, from a QCD sum 
rule point of view, it is much easier to form an state, separated from the continuum, with three quarks 
than with four quarks. 
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